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Abstract. We analyze monotone finite difference schemes for strongly degen- 
erate convection-diffusion equations in one spatial dimension. These nonlinear 
equations are well-posed within a class of (discontinuous) entropy solutions. 
We prove that the L 1 error between the approximate and exact solutions is 
©(Ax 1 / 3 ), where Aa; is the spatial grid parameter. This result should be 
compared with the classical 0(Ax 1 ^ 2 ) error estimate for conservation laws 
1221 , and a recent estimate of ©(Ax 1 / 11 ) for degenerate convection-diffusion 
equations 1201 . 
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1. Introduction 

Nonlinear convection-dominated flow problems arise in a range of applications, 
such as fluid dynamics, meteorology, transport of oil and gas in porous media, 
electro-magnetism, as well as in many other applications. As a consequence it has 
become a very important undertaking to construct robust, accurate, and efficient 
methods for the numerical approximation of such problems. Over the years a large 
number of stable (convergent) numerical methods have been developed for linear 
and nonlinear convection-diffusion equations in which the "diffusion part" is small, 
or even vanishing, relative to the "convection part" of the equation. There is a 
large literature on this topic, and we will provide a few relevant references later. 

One central but exceedingly difficult issue relating to numerical methods for 
convection-diffusion equations, is the derivation of (a priori) error estimates that 
arc robust in the singular limit as the diffusion coefficient vanishes, avoiding the 
exponential growth of error constants. This problem has been resolved only partly 
in special situations, such as for linear equations or in the completely degenerate 
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case of no diffusion (scalar conservation laws). For general nonlinear equations 
containing both convection and (degenerate) diffusion terms this is a long standing 
open problem in numerical analysis. 

This paper makes a small contribution to this general problem by deriving an 
error estimate for a class of simple difference schemes for nonlinear and strongly 
degenerate convection-diffusion problems of the form 

d t u + d x f(u) = d 2 x A(u), (x,t)en T , 

u(0,x) = u°{x), x G R, 

where Ily =1x (0, T) for some fixed final time T > 0, and u(x,t) is the scalar 
unknown function that is sought. The initial function uq : R — > K is a given 
integrable and bounded function, while the convection flux / : K — s- R and the 
diffusion function are given functions satisfying 

/, A locally C ; A(0) = 0; A nondecreasing. 

The moniker strongly degenerate means that we allow A'(u) = for all u in some 
interval [a, (3] C K. Thus, the class of equations becomes very general, including 
purely hyperbolic equations (scalar conservation laws) 

d t u + d x f(u) = 0, (1.2) 

as well as nondegenerate (uniformly parabolic) equations, such as the heat equation 
dtu = d x u, and point-degenerate diffusion equations, such as the heat equation with 
a power-law nonlinearity: d t u = d x (u m d x u), which is degenerate at u = 0. 

Whenever the problem is uniformly parabolic (i.e., A' > rj for some n > 0), 
it is well known that the problem admits a unique classical (smooth) solution. 
On the other hand, in the strongly degenerate case, (|f ,f [) must be interpreted 
in the weak sense to account for possibly discontinuous (shock wave) solutions. 
Regarding weak solutions, it turns out that one needs an additional condition, 
the so-called entropy condition, to ensure that (jTTTJ) is well-posed. More precisely, 
the following is known: For uq G L 1 (K.) n L°°(R), there exists a unique solution 
u G C((0,T);L 1 (R d )), u G L°°(II T ) of (TTl]) such that d x A(u) G L 2 (T1 T ) and for all 
convex functions S : R ->• R with q' s = f'S' and r' s = A'S', 

d t S(u) + d x q s {u) - d 2 x r s {u) < in the weak sense on [0, T) x R. (f .3) 

The satisfaction of these inequalities for all convex S is the entropy condition, 
and a weak solution satisfying the entropy condition is called an entropy solution. 
The well-posedness of entropy solutions is a famous result due to Kruzkov [3T| 
for conservation laws (jf .2[) . and a more recent work by Carrillo [5] extends this 
to degenerate parabolic equations (jf .f [) . These results are available in the multi- 
dimensional context, and we refer to [UHHI f° r & n overview of the relevant literature. 
For uniqueness of entropy solutions in the BV class, see [261 128) . 

One traditional way of constructing entropy solutions is by the vanishing viscosity 
method, which starts off from classical solutions to the nondegenerate equation 

dtu v + d x f{u n ) = dlA(u v ) + r]dlu v , n > 0, 

and establishes the strong convergence of as ry — > by deriving BV estimates 
that are independent of 77, see Vol'pert and Hudjaev [27] . 

Besides proving that u v converges in the L 1 norm to the unique entropy solution 
u of (jl.ip . it is possible to prove the error estimate 

\\u v {-,t) -u(-,t)\\ Ll < Cy/rj, (whenever u$ G BV), (1.4) 

see [T3] (cf. also [H]). The error bound (|1.4[) can also be obtained as a consequence 
of the more general continuous dependence estimate derived in [5] , see also [SI [H] • 
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Herein we are interested in the much more difficult problem of deriving error 
estimates for numerical approximations of entropy solutions to convection-diffusion 
equations. Convergence results (without error estimates) have been obtained for 
finite difference schemes [12] (see also [Ml [19]); finite volume schemes [16] (see also 
[5]); operator splitting methods [T7]; and BGK approximations pOH], to mention 
just a few references. For a posteriori estimates for finite volume schemes, see [23] . 

To be concrete in what follows, let us for simplicity assume /' > and consider 
the semi-discrete difference scheme 

where Uj(t) w u(t,jA.x) and Ax > is the spatial mesh size. Convergence of 
this scheme can be proved as in the works [121 113] , where explicit and implicit time 
discretizations are treated. Denote by mai(i, t) the piecewise constant interpolation 
of {uj(t)}j. The basic question we address in this paper is the following one: Does 
there exist a number r £ (0, 1) and a constant C, independent of Ax, such that 

\\u Ax (;t)-u(;t)\\ Ll <CAx r , (1.6) 

where u is the unique entropy solution of (jl.ll) . We refer to the number r as the 
rate of convergence. 

In the purely hyperbolic case (|1.2[) (A' = 0), the answer to this question is a 
classical result due to Kuznetsov [22], who proved that the rate of convergence is 
1/2 for viscous approximations as well as monotone difference schemes, and this 
is optimal for discontinuous solutions. The work of Kuznetsov [22] turned out to 
be extremely influential, and by now a large number of related works have been 
devoted to error estimation theory for conservation laws. We refer to [7] for an 
overview of the relevant results and literature. 

Unfortunately, the situation is unclear in the degenerate parabolic case (|1.1[) . Let 
us expose some reasons why adding a nonlinear diffusion term to (|1.2j) can make the 
error analysis significantly more difficult than in the streamlined Kuznetsov theory 
First of all, it is well known that the purely hyperbolic difference scheme 



has as a model equation the second order viscous equation 

d t u + d x f{u) = ^-d 2 j{u), 

an equation that is compatible with the notion of entropy solution for (|1.2|) . Indeed, 
an error estimate for this viscous equation is highly suggestive for what to expect 
for the upwind scheme (|1.7p (this is of course what Kuznetsov proved). However, 
for convection-diffusion equations such as (|1.1[) the situation changes. The model 
equation for (|1.5[) is no longer second order but rather fourth order: 

d t u + d x f(u) = d 2 x A(u) + ^d 2 J{u) - ^-d 4 x A(u); 

hence the error estimate (jl.4l) appears no longer so relevant for numerical schemes. 
Another added difficulty comes from the necessity to work with an explicit form 
of the parabolic dissipation term associated with p. II) . Indeed, in the analysis one 
needs to replace (|1.3[) by the following more precise entropy equation [5] 

d t \u-c\+d x (sign(u-c)(f(u)-f(c))-dl\A(u)-A(c)\ 

9 (1-8) 
= -sign'(A(u) - A(cj) \d x A{u)\ 2 , c G K, 

which is formally obtained multiplying (|1.1[) by sign (A(u) — A(c)), assuming for 
the sake of this discussion that A'(-) > 0. The term on the right-hand side is 
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the parabolic dissipation term, which is a finite (signed) measure and thus very 
singular. To illustrate why the parabolic dissipation term is needed, let u(y,s) and 
v(x,t) be two solutions satisfying (|1.8[) . In the entropy equation for u(y,s) one 
takes c = v(x,t), while in the entropy equation for v(x,t) one takes c = u(y,s). 
Adding the two resulting equations yields 

(fit + d s ) \u-v\ + (d x + a y )(sign(u - v)(f(u) - /(«)) 

- {dl + d 2 y ) \A(u) - A(v)\ = -sign'(A(u) - A(v))(\d y A(u)\ 2 + \d x A(v)\ 2 ), 

By adding —2d^. y \ A(u) — A(v)\ to both sides of this equation, noting that 

-2d 2 xy \A{u) - A(v)\ = 2sign'(A(u) - A(v))d y A(u)d x A(v), 

we arrive at 

(dt + d s ) \u-v\ + (d x + a,)(sign( U - v)(f(u) - f(v)) 
-(dl-2dl y + dl)\A{u)-A{v)\ 

(19) 

- -sign' (A(u) - A(v))(\d y A(u)\ - \d x A(v)\f 
<0, 

from which the contraction property -4 \\u(-, t) — v(-, t)\\ L1 < follows [5]. Similarly, 
to obtain error estimates for numerical methods, it is necessary to derive a "discrete" 
version of (|1.9[) with v replaced by ua x - The main challenge is to suitably replicate 
at the discrete level the delicate balance between the two terms in (jl.9|) involving 
A; the difficulty stems from the lack of a chain rule for finite differences. 

Despite the mentioned difficulties, we will in this paper prove that there exists 
a constant C, independent of Ace, such that for any t > 0, 

||ua*(->*)-«(->*)IIli ^ CAxi. 

The only other work we are aware of that provides L 1 error estimates for numerical 
approximations of (jl.ip is ^U\; therein (|1.6p is established with r = jj; if A is 
a linear function, then the convergence rate is the usual one, namely r = |. In 
addition to the semi-discrete scheme (jl.7[) . we will prove similar results for fully 
discrete (implicit and explicit) difference schemes. 

Roughly speaking, the reason is two-fold for why we can significantly improve 
the result in [50]. First, we are herein able to provide a more faithful analog of 
(|1.9[) at the discrete level. Second, since sign'(-) is singular, one has to work with 
a Lipschitz continuous approximation sign E (-) of the sign function sign(-). The use 
of this approximation breaks the symmetry of the corresponding entropy fluxes, 
and introduces new error terms that depend on the parameter e; the process of 
"balancing" terms involving Ax and e lowers the convergence rate (to r = jj) [50] . 
In the present paper we are able to dispense with this balancing act. Indeed, we 
show that it is possible to send e — > independently of Ace. 

The remaining part of this paper is organized as follows: In Section [2] we list 
some relevant a priori estimates satisfied by viscous approximations and entropy 
solutions, and provide a definition of entropy solutions. The semi-discrete difference 
scheme is defined and proved to be well-posed in Section [31 We also list several 
relevant a priori estimates. Section 2] is devoted to the proof of the error estimate. 
In Section [S] we show that the proof in Section 0] can be adapted to a fully discrete 
scheme that is implicit in the time variable. In fact, we go through all the steps of 
the proof and provide the details where there are considerable differences between 
the two cases. In Section [S] the explicit version of the scheme is treated. We end 
the paper with a few concluding remarks in Section [7] 
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2. Preliminary material 

Set := A(u) + rju for any fixed r\ > 0, and consider the uniformly parabolic 

problem 

f«? + /(«")« = A"(«") XX> (x,i)en T , 

1 ^(x.O) = u°(aj), xeR. 

It is well known that (|2.1I) admits a unique classical (smooth) solution. 

We collect some relevant (standard) a priori estimates in the next three lemmas. 

Lemma 2.1. Suppose u° G L°°(R) n L X (R) n W(R) 7 and Zet u 17 oe the unique 
classical solution of (|2.1[) . TTien /or any i > 0, 

^''(''^IIl^R) - II U °IIl1(R) ' 

IK(->*)lli°°(R) < II«°IU«(r), 

I u,? ("j*)Ibv(k) - I u °Ibv(r)- 

For a proof of the previous and next lemmas, see for example |27| . 

Lemma 2.2. Suppose u° G L°°(R) n L X (R) and /(it ) - ^(u )^ G BV(R). Let vP 
be the unique classical solution of (|2.1[) . Then for any t\,t2 > 0, 

||«"(.,* 2 ) - u^^h)^ < \f(u°) ~ A(u°) x \ Bvm \t 2 
Regarding the following lemma, see [UJ H2] ■ 

Lemma 2.3. Suppose u° G J L 00 (R)nL 1 (R) and /(it ) - A(u°) x E L°°(R)nW(R). 
Let it'' be the unique classical solution of (|2.1j) . T/ien /or any t > 0, 

||/(«*M))-A(u*M))j LOO(R) < \\f(u )-A(u°) x \\ LooW , 

\f(u^.,t))-A(uV(.,t)) x \ Bvm < \ f(u°)-A(u°) x \ Bvm . 

Note that \\A(u n ) x \\ Lal ^ LCX ,^ and || A(it ?? ) 2 ; :E |j L=0 ^ i i^ are bounded independently of 
77 provided that A(u°) x is in W(R). 

The results above imply that {u v } v>0 is relatively compact in C([0, T]; L 1 1 oc (R)). 
If u = lim^o it'', then 

IK'-«ll L1( n r) <Cv 1/2 , 
for some constant C which does not depend on 77, see fTJ]. Moreover, it is an entropy 
solution according to the following definition. 

Definition 2.1. An entropy solution of the Cauchy problem is a measurable 
function u = u(x, t) satisfying: 

(D.l) ue L oo (n T )nC((0,r);L 1 (K)). 

(D.2) A(u) G ^((OjTJjlT^R)). 

(D.3) For all constants c G R and test functions < ip G C^°(R x [0,T)), the 
following entropy inequality holds: 

|it - c| v?t + sign (it - c) (/(«) - f(c))tp x + \A(u) - A(c)| deedt 

\uq — c\ <p(x, 0) dx > 0. 



The uniqueness of entropy solutions follows from the work [5] . Actually, in view 
of the above a priori estimates, the relevant functional class is BV(Ut), in which 
case we can replace (D.2) by the condition A(u) x £ L°°(Ht). For a uniqueness 
result in the BV class, see 1251. 
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3. Difference scheme 

We start by specifying the numerical flux to be used in the difference scheme. 

Definition 3.1. (Numerical flux) We call a function F e C 1 (IR 2 ) a two-point 
numerical flux for / if F(u, u) = f(u) for net. If 

d d 
— F(u,v)>0 and ^F(u,v)<0 
ou ov 

holds for all a,uf!l, we call F monotone. 

Let F u and F v denote the partial derivatives of F with respect to the first and 
second variable, respectively. We will also assume that F is Lipschitz continuous. 
Let Ax > and set Xj = jAx for j € Z, and define 

D ± a j =± <Tj±1 ~ <Tj 



Ax ' 
for any sequence 

We may now define a semi-discrete approximation of the solution to (|I.1[) as the 
solution to the (infinite) system of ordinary differential equations 

j&u^+D.Fj+y^D-D+Aiuj), t>0, 

\«j(o) = ^4«°(*)«fa. 3 [ j 

where Fj +1 / 2 = F(uj, Uj+i) is a numerical flux function and Ij = (£j_i/2, 2^+1/2]- 
The problem above can be viewed as an ordinary differential equation in the 
Banach space ^ 1 (Z) (see, e.g., [131 )■ To get bounds independent of Ax we define 

Hi = Ax ^2\ a j\ and W\bv = l^ 1 ~ = H^+^lli • 
j 3 

If these are bounded we say that a = is in i 1 and of bounded variation. 

Let u(t) = {uj(t)}. &z , u° = {uj(0)}j ez , and define the operator A : i 1 — > i 1 by 
(A(u))j := D-(F(v,j,Uj+i) - D + A(ujj). Then (j3"Tj) takes the following form 

^+A(u) = 0, t>0, u(0)=u°. 
dt 

This problem has a unique continuously diffcrcntiablc solution since A is Lipschitz 
continuous for each fixed Ax > 0. This solution defines a strongly continuous 
semigroup S(t) on i 1 . If S also satisfies 

WStyu-StyvW-L < for u,vef, 

we say that it is nonexpansive. The next lemma sums up some important properties 
of the solutions to (|3.ip (for a proof see |llj). 

Lemma 3.1. Suppose that F is monotone. Then there exists a unique solution 
u = {uj} to (|3.Ij) on [0, T] with the following properties: 

(a) ||«(t)||i < \\u%. 

(b) For every j G Z and t G [0, T], 

infK} <«,•(*)< sup K}. 

k k 

(c) |«(t)| BV < |«°| BV . 

(d) If v = {vj} is a another solution with initial data v° then 

H«(t) -«c*)lli < ||«° - « IU - 
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Lemma 3.2. If F is monotone, then 

\\F( Uj ,u J+1 ) - D + A( Uj )\\^ < - D + A(u°)\\ eoo , (3.2) 

\F( Uj , u j+1 ) - D + A( Uj )\ BV < \F(u%u° ]+l )-D + A{u a )\ BV . (3.3) 

Furthermore, t > {uj(t)} jgZ is t Lipschitz continuous. 

Proof. The proof follows [T2]. Let Vj = AxJ2k<j ^lt- Then Vj satisfies 
o 

Vj=Ax D-(D + A(u k )-F(u k ,u k+1 ))=D + A(u j )-F(u j ,u j+1 ), (3.4) 



and we may define Vj for all t £ [0,T]. Note that {vj(t)} is in I 1 for all t by 
Lemma I3~T1 Differentiating p.4[) with respect to t we obtain 



cfoj 1 
dt ~ Ax 



du j+1 duj 

a ( u 3+i)— — a ^'~dT 

duj „ . . du 



eft v ■" dt 



F u (uj,u j+1 ) — - F v (uj,u j+ i)- 



dt dt 
where a(w) = A'(u). Note that D-Uj = and -D+Vj = d "j t +1 . Therefore 



^-a(uj) + F u (wj,w i+ i)^ (3.5) 

Assume Vj (to) is a local maximum in j. Then D + Vj a (to) < and D-Vj (to) > so 
-4^-(to) < since i* 1 is monotone. Similarly, if Uj (to) is a local minimum in j, then 
^(t ) > 0. Then inequality flUD) follows by the fact that {u,(t)} € ^. Consider 
([3TB]1 . We want to show that ^ (KOIw) < 0. Now, 

3 3 



so we may use (|3.5[) . Thus 

|k*)Ibv 

= XI ( -K- a ( u 3+z) ~ F v(uj+i,u j+2 ) J (D+w J+ i)sign(w J+1 - Vj) 



'+ ^ I 



Ax 

(^ a (%"+i) +i;l «( w i+i' w J+2)J |-Dh 
+ (^ a( ^ + ((-D-«i)sign(u,-+i - Vj)) 



j 

<0, 



8 K. H. KARLSEN, N. H. RISEBRO, AND E. B. STORR0STEN 

since a(u) > 0, F v < 0, and F u > 0. Given the preceding estimates, the I 1 Lipschitz 
continuity is straightforward to prove. □ 

It turns out that we need more conditions on F than mere monotonicity. 

Definition 3.2. Given an entropy pair (ijj, q) and a numerical flux F, we dchnc 
Q £ C 1 ^ 2 ) by 

Q(u, u) = q(u), 

d d 
—Q(v,w) = ip (v)—F(v,w), 

d d 
—Q(v,w) = ip'(w) — F(v,w). 
aw ow 

We call Q a numerical entropy flux. 

The next lemma gives a sufficient condition on the numerical flux to ensure that 
there exists a numerical entropy flux. 

Lemma 3.3. Given a two-point numerical flux F, assume that there exist C 1 
functions F\ , F2 such that 

F(u,v) = F 1 (u)+F 2 (v), F[{u)+F^u) = /'(«), (3.6) 

for all relevant u and v. Then there exists a numerical entropy flux Q for any 
entropy flux pair [ip, q). 

Proof. Let (tp, q) be an entropy pair. Then q has the form 

1/(z)f'(z)dz + C, 

for some constant C. Define Q by 

tp'(z)F[(z)dz + J il}'{z)F! 2 {z)dz + C. (3.7) 
It is easily verified that Q is a numerical entropy flux. □ 

Let us list a few numerical flux functions to which Lemma 13.31 applies. 
Example 3.1 (Engquist-Osher flux). Let 

f + (s) = max(/'( S ),0) and /1(a) = min(/'(a), 0). 
Then, in the terminology of Lemma |3~31 let F(u,v) = Fi(u) + F2(v) with 

F 1 (u) = f(0)+[ f'(s)ds and F 2 (v) = [ f'_ (s) ds. 



It is easily seen that the criteria given in Lemma 13.31 are satisfied, and F is also 
clearly monotone. 

Example 3.2. Let a, b £ R and define 

= af{u) + bu and F 2 (v) = (I - a)f(v) - bv. 
Note that F(u,v) = F±(u) + F 2 (v) is monotone if 

a inf {/'(«)} > -b and (I - a) sup{/'(u)} < b. 

u u 

This example includes both the upwind scheme and the Lax-Friedrichs scheme. 
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From a more general point of view we may consider any flux splitting, that is, 
f(u) = f+(u) + f-(u) with (f + )'(u) > and (/")'(«) < for all ueE. Then the 
numerical flux 

F(u, «) = /+(«) +f-(v) 
satisfies the assumptions of Lemma 13.31 Note also that any convex combination of 
numerical flux functions which satisfy the hypothesis of Lemma 13.31 itself satisfies 
the assumptions of the lemma. 

If p.6[) holds, then we have a representation of Q given by (|3.7p . It follows that 

Q(u,v) = q(u)+ [ ip'(z)F^(z)dz. (3.8) 



Note that we may obtain another representation depending on F\ by splitting up 
the first integral. 

4. Error estimate 

Let {%}j £Z be the solution to (|3.1[) . We associate with it the piecewise constant 
function 

UAx(x,t) = Uj(t) for x G Ij. (4-1) 
To derive the error estimate we need many of the uniform bounds from Sections [2] 
and [3l For these estimates to hold independently of Ax, we make the following 
assumptions on the initial data u . 

(i) u° G L 1 (R) n L°°(R) n BV(R). 

(ii) A(u°) x G BV(M). 

We may now state the theorem. 

Theorem 4.1. Let u be the entropy solution to (jl.ip and {uj(t)}. gZ so/we the 
semi-discrete difference scheme (|3.ip . I/it satisfies (i) and (ii) above, then for all 
sufficiently small Ax, 

||«A«('>*)-«('>*)ll£i(R) < |hL-«°|| L i (R) +^TAx3, ig [0,T], 

where the constant C't depends on A, /, u , and T, but not on Ax. 

Let us define some of the functions we are going to work with. First, we will use 
the following approximation of the sign function: 

. fsin(lf) for \a\ < e, 

sign, {(J) = i . f . 

I sign (a) otnerwise, 

where e > 0. Note that sign 6 is continuously differentiable and non-decreasing. We 
define 

\u\ e = / sign e (z) dz. 
Jo 

Furthermore, we introduce an entropy pair (ip e , q e ) defined by 



ip e (u, c) = J sign £ (A(z) - A(c)) dz, 

pU pU 

q E (u,c)= ^' E {z,c)f{z)dz= sign s (A(z)-A(c))f'(z)dz, 



where is the derivative with respect to the first variable. 

Lemma 4.1. Suppose A' > 0. Let u = u(y,s) be the classical solution of (|l.ip . 
Then for any constant c 6 R, 

d s 4> e (u,c) + d y q e (u,c) - d%\A(u) - A(c)\ £ = -d y ip' E {u,c)d y A(u). 
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Proof. Multiply equation ([l.ip by ip' £ (u,c) to obtain 

d s ip e (u,c) + d y q E {u, c) = ip' E (u, c)dyA(u). 
The term on the right may be rewritten according to 

dyMfa c)d y A(u)) = d y ^' e (u, c)d y A(u) + i>' e {u, c)d 2 y A{u). 
By the chain rule 

dM{u,c)d v A{u))=dl\A{u)-A{c)\ e . 
Combining these equalities proves the lemma. □ 

The next lemma is a simple identity taken from [8] . 
Lemma 4.2. For any differ entiable function g and all real numbers a,b,c, 

^{a,c){g(b)-g(a)) = [ ^(z,c)g'(z)dz - f " ^' s (z,c)g'(z)dz 



+ [ ^(z,c)(g(z)-g(b))dz. 

J a 

Proof. Integration by parts yields 

#(C, c)(g(C) - g(b)) = J ifaz, c)g'(z) dz + J ^(z, c)(g{z) - g(bj) dz, 

for any ( e K. Take the two equations obtained by taking £ = a and C = b and 
subtract one from the other. □ 

Lemma 4.3. Let uj be the solution to (|3.1[) . Then for all c el 

d t ip e (uj,c) + D^Q c (uj,u j+ i) -D-D + \A( Uj ) - A(c)\ e 
1 



< -- 



(Ax) 2 
1 



^(z,c)(A(z) - A{u j+1 )) dz 

u j + 1 



^{z,c){A{z)-A{u j - l ))dz, 



(Ax) 2 

where Q c (u, v) := Q\(u) + Q 2 (v), 

Qt(u) := J ^ s (z,c)Fi(z)dz, Q c 2 (v) := J tf s {z,c)F* 2 (z)dz, 

for all real numbers u and v. 
Proof. From (|3.1|) it follows that 

il)' e (uj,c)dtUj + ip' E (uj,c)D_F(uj,Uj + i) = iJ/ £ (uj,c)D-D + A(uj). 
Note that 

^'e{uj,c)D_F(uj,u j+1 ) = ij' s (u j ,c)D^F 1 (u j ) + i(j' E (u j ,c)D + F 2 (u j ), 
and so we may apply Lemma 14.21 Let g = F\ . Then we obtain 

^(t^cJD.JKu,-) = D_Qi( Uj ) 1 ^(z,c)(F 1 (z)-F 1 (u j . 1 ))dz. 

Similarily, let g = F 2 to obtain 

^ e (nj,c)D + F a (uj) = D + Q c 2 (u 3 ) + ^ £ ^(z,c)(F 2 (z) - F 2 (u J+1 ))dz. 

Finally, apply lemma l4~2l twice with g = A. Adding the equations we obtain 
i)' e ( Uj , c) (A(u^i) - 2A(uj) + A(u j+1 )) 
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ip' e (z,c)A'(z)dz+ / ip' E (z,c)A'(z)dz 



+ / ^(z,c)(A(z)-A(u j+1 ))dz 



+ i//J{z,c){A(z)-A(u j - 1 ))dz. 



Note that 



%j)' £ (z,c)A'(z)dz + \ ip' e (z,c)A'(z)dz 



9 



_|A(z)-A(c)| £ dz + — |A(z)-A(c)| e <fe 



9z 



= [|^K-i) - A{c)\ e - 21^-) - A(c)| e + \A(u j+1 ) A(c)\ e ] . 
Combining the above computations we obtain 

d t ip e (uj,c) + D_Q c (uj,u j+ i) - D-D + \A(uj) - A(u)\ e = -E c (uj-x,Uj,Uj +1 ), 
where 

1 



E c (uj-i,Uj,u j+ i) 



^'(z,c)(Ji(«) - Fxiuj-^dz 



l r 3 

— ^(z,c)(F 2 (z)-F 2 (u j+1 ))dz 

^ x Ju j + 1 



1 



(Ax) 2 J Uj+1 

1 



^(z,c)(A(z)-A(u j+1 ))dz 



(Ax) 2 



^(z,c)(A(z)-A(u j . 1 ))dz. 



The result follows from the monotonicity of F. 



□ 



We shall need the next lemma, which deals with a mixed term, in order to carry 
out the "second order" doubling-of-the- variables argument. 

Lemma 4.4. Let {uj} be some sequence and u some differentiable function of y. 
Then 



Ax 



ip' £ (z, u) dz + 



i/)"(z,u) dz J d y A{u) 



-d v (D_ + D+)\A(uj) - A{u)\ £ 



Proof. Let a, & be hxed real numbers. Then 

" ^(z,u)dzA(u) v 



signl(A(z) - A(u))A(u) y A'(z) dz 



0_ 

dy 

d_ 

dy 



■— I J sign E (A(z)-A(u))A'(z)dz 



(\A(b)-A(u)\ e -\A(a)-A(u)\ B ). 



Let a = Uj-i, b = Uj and a = Uj, b = Uj+i. Then add up the resulting equations 
and divide by Ax. □ 
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We are now in a position to carry out the doubling-of-thc- variables argument. 

Lemma 4.5. Suppose A' > 0. Let u = u(y, s) be the classical solution to and 
let {uj} = {uj(t)} be the solution to (|3.1[) . Then 

d t il) s {uj t u) + d s ip £ (u, Uj) + d y q e (u,Uj) + D_Q u (uj,Uj + i) 

- (d 2 y + d v {D_ + D+) + D_D + )\A{u 3 ) - A(u)\ £ ) < -Ej, 

where := E E [u](uj-i, uj, Uj+i) with 

1 f U3 

E e [u]{uj-i,Uj,u j+1 ) := I i> e {z,u){A(z) - A{u j+l ))dz 



u j+i 



ip"{z, u)(A(z) - A(uj-i)) dz 



"j-i 



(Ax) 
1 [ U] 

— I ip"(z,u) dzd y A{u) 
+ 1 

■0" (z, u) dz d y A(u) 



Ax 

+ d v tp' E (u,Uj)d y A(u). 

Proof. Let c = Uj in Lemma 14.11 and c = u in Lemma 14.31 Then add up the 
equations together with Lemma T4. 41 □ 

Remark 4.1. Note that E^ is a function of y, s, t. 

In what follows it will be necessary to work with the piecewice constant approx- 
imation dchncd in (|4.1j) . To do this we introduce some new notation. Let the shift 
operator 5*0- be defined for any if : TLt — > R by 

S a ip(x,t) = ip(x + a,t), 

and the difference quotient be defined by 

S±AxP - V 

D±lf = Ax 

Note that for any two functions u,v of x we have D+(uv) = SaxuD+v + (D + u)v. 
If uv has compact support it follows that 



/ (D+u)vdx = — / uD-vdx. 
Jr Jr 



We will refer to these identities as the Leibniz rule for difference quotients and 
integration by parts for difference quotients. We will frequently integrate over the 
domain 11^ . To avoid writing four integral signs we will in general write one for 
each domain ILp and let dX = dxdtdyds. 

Lemma 4.6. Suppose A' > 0. Let u^x — UAx(x,t) be defined by (|4.1[) . and let 

u = u(y,s) be the classical solution of . Let p € Co°(R) satisfy 

supp(p) C [-1, 1], p(-o-) = p{a), p{a) > 0, J p{a) do = 1, 

and set 

W r (x) = -p(-), P a (£) = -p(-), Pr (t) = — P ( — ) , 

r \r J a \a J fo\ r o/ 

for positive (small) r, a and r$. Let v and r be such that < v < r < T and define 

^ a {t):=H a (t-v)-H a (t- T ), H a (t)= [ Pa (0dt 
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Let 

ip(x,t,y,s) = ip a (t)uj r (x - y)p ra (t - s). 
To ensure if\t=o = 0, if\ s= o = 0, we choose v and r such that < ro < min(v, T — t) 
and < a < mm.{y — ro, T — r — ro)- Then 



IL 



\UAX - u\ p a (t - v)uJ r Pr dX 

+ J J ^ sign (u Ax - u) (f(u Ax ) - f(u)) (D + ip + tp y ) dX 
+ // / sign (z-u)F^(z)dz ) D+tpdX (4.2) 

+ JJ 2 \A(uax) - M u )\ ( D - D +V + ( D + + D -)<Py + Vyy) dX 



> 



1 1 \uax - u\ p a (t - T)oj r p ra dX + lim mf / / E Ax tpdX, 



where E% x (x, t, y, s) = E £ 3 (t, y, s) for x G Ij. 

Remark 4.2. Note that both 

<Px + <Py — and ip xx + 2ip xy + tp yy = 0. 

In equation (|4.2[) these expressions appear with difference quotients instead of x- 
derivatives. We expect that these equalities turns into good approximations as long 
as Ax tends relatively fast to zero compared to r. We will show that this is the 
case in what follows. 

Proof. By Lemma T4. 5 1 it follows that 

9 t tpe(UAx,u) +d s 1p £ (u, UA X ) + d y q £ {u,UA x ) + D-Q u (uAxiSaxUAx) 

- (dl + d y (D_ + D+) + D_D + )\A(uax) - A(u)\ e < -E% x , 

for all (x, t, y, s) G Ily. Let us multiply with ip and integrate over n T . Using both 
ordinary integration by parts and integration by parts for difference quotients, we 
obtain 

1pe{u*Ax,u)(pt + 1p e (u, U A x)V>s dX 



+ J J 2 qe(u,UAx)fy + Q U (uAx,SAxUAx)D + tfdX 

+ JJ^ \A(uAx)-A(u)\ s (ip yy + {D_+D + ) Vy + D_D + ip)dX 
> JJ n2 E Ax <PdX, 



We want to take the limit as e I 0. Consider the first term on the left. By the 
dominated convergence theorem, for any a, b € K, 

linage (a, b) = lim / sign £ (A(z) - A(b)) dz = \a - b\ , 

since A' > 0. It follows that 

lim ip e (u Ax, u) = lim ip E (it, uax) = |«Ax — «|- 

£10 elO 

Furthermore, 

((ft +(p s )(x,t,y,s) = (p a (t - v) - p a (t - r))uj r (x - y)p ro (t- s), 
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so by the dominated convergence theorem 



lim// ip £ (uA x ,u) ( Pt+'<Pe(u,UA X )<PsdX 

EiO J J U 2, 

\u Ax - u\ p Q (t - v)u) r p ro dX - \u Ax - u\ p a (t - r)uj r p ro dX 

JJn 2 T 

Consider the second term on the left. By (|3.8[) we obtain 

Q v (uax, SaxUAx) = q e (u A x, it) + / sign e (A(z) - A^l^F^z) dz. 

J MAa: 

Since A 1 > 0, 

limq e (uAx,u) = lim / sign £ (A(z) - A(u Ax ))f'(z) dz 
= sign (u Ax - u) (f(v,Ax) - f(u))- 



It follows that 



As above 



\imQ u {u Ax ,S Ax u Ax ) = sign(w A 2- - u) {f{u Ax ) - f(u)) 

+ / sign (z - u) F^z) dz. 



limq e (it, u Ax ) = sign (u - u Ax ) (f(u) - f{u Ax )) 

ei.0 

= sign (uax - u) {f{u Ax ) - f(u)). 
Hence, again by the dominated convergence theorem, 

lim // q e (u, u Ax )<p y + Q u (u Ax , S Ax it Ax )D + <p dX 
^° J Jul 

sign (u Ax - u) (f(uAx) - f(u))(<p y + D + ip) dX 



/ Slf 



' fSAi"AT \ 

/ sign (z — u) F' 2 (z) dz D + ip dX. 



Lemma 4.7. Let E Ax and tp be defined in Lemma \J h .b\ Then 

liminf // E Ax ipdX > / liminf I / E Ax ipdyds\ dxdt. 



Proof. Let 
and 



f s (x,t) := / E e Ax ipdyds 



□ 



h £ (x,t):= / d y (D- +D + )\A(u Ax ) - A(u)\ e <pdyds. 

Recall that E Ax (x,t,y,s) = Ej(t,y,s) for x <G Ij, where Ej is defined in Lemma l4~5l 
Note that 



E$>--£- I ' i>' e \z,u)dzd y A{u) 
_ _L P +1 ^( z ,u)dzd y A(u), 



ERROR ESTIMATES FOR CONVECTION-DIFFUSION EQUATIONS 



IS 



so by Lemma Eni it follows that f £ > h £ . Using integration by parts and the triangle 
inequality we obtain the bound 

\h e \ < (\D_A(u Ax )\ + \D+A{u Ax )\) Qf \ip y \ dyds^j =: h. 

It follows by Lemma 13.11 that h is an integrable nonnegative function such that 
— h < f e - By Fatou's lemma we obtain 



liminf / f^dxdt> / liminf f s dxdt. 
40 Jur Jur 40 



□ 



Note that as e i the terms in Ej concentrate on the domains specified by 
u G int (it,-, itj-|_i), u G int(wj_i, Uj), or u = uj. In order to analyze this limit we 
will need the following elementary lemma: 

Lemma 4.8. Let {uj}j £Z be some sequence in K and let A : K — > M a strictly 
increasing continuously differ entiable function. For any u G K t/iere exisi sequences 
{ r j t }- gZ i {^Ijgz s«c/i i/iai /or eac/i j G Z ooi/i and 0^ are m int(u, ■, Uj±i) 
and 

D±sign e (^(«j) - i4(u)) = sig n ;(A(r±) - A(u))£)±A(^), 
D±|A( Ui ) - A(u)\ £ = sign £ (A(9f) - A(u))D ± A{ Uj ). 
If u is a differentiable function of y then for each j G Z, 

sig<(A(r±) - - -(sign e (A(0±) - A(u))) y . (4.3) 

Bo£/i {tJ^}^^ and . gZ depend on u and e. 

Proof. The first statement is a direct consequence of the mean value theorem. 
Consider (|4.3[) . First note that t~ = rjt_i and 0J = OjLi, so it suffices to consider 
and If uj = uj + i then 8j = tj is independent of u and hence of y, so (|4.3p 
follows by the chain rule. In general, 

sigri £ (A(Tj) - A(u))A(u) y D + A(uj) = D + sign e (A(uj) - A(u))A(u) y 

= -D+(\A( Uj ) - A(u)\ e ) y 

= -sign £ {A(6 3 ) - A(u)) y D + A(u 3 ). 

In the case Uj 7^ %+i we have D + A(u,j) ^ and (|4.3[) follows. □ 

The following result is concerned with the pointwise limit of sign e (A(9j) — A(u)) 
as £ | 0. The explicit formula for this limit, which will be used later, shows that the 
limit is in fact a Lipschitz continuous function in the case that A(uj) ^ A{uj±\). 

Lemma 4.9. Let 

{la— crl — \b— a - / i 

sign(a -a) if a = b, a ^ b. 
if a = b = a, 

for any real numbers a and b. Under the same assumptions as in Lemma 
1 ™sign e (A(6»= t ) - A(u)) = -sg {A[U])tA{U]±l}} (A(u)) . 

Furthermore, if a ^ b then 

-1 if a < min {a, b}, 

Sg( a ,6) O) = { T^Ff (<r~ U b + a )) l f ae int(a,&), 

1 if a > max {a, b}. 
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Proof. To prove the first statement we consider the case of 9~j~ . The same argument 
applies to 9 J . Recall the definition of 9^ : 

sign, (A(9+) A(u)) (A(u j+1 ) A( Uj )) = \A( Uj+1 ) A{u)\ e - \A( Uj ) A(u)\ e . 

If Uj + i = Uj, then 9j~ = Uj for all u and e, since 9^ 6 int(wj, Wj+i). Thus in this 
case 

]imsignJA(ei)-A(u)) = 1° if « = «,-, 

<40 " J I sign (^(lij) — otherwise. 

Now assume that ^ 0. Then 

\A{u ]+1 )-A{u)\ e -\A(u 3 )-A{u)\ £ 



sign E (A(9+) - A(u)) - 

3 Mui+i) - A(uj) 



and the result follows by letting e I 0. Let us prove the second statement. First 
observe that all expressions are symmetric in a and 6, so we may assume that a < b. 
Under this assumption we have 

- a)sg( a ,&) 0) = \a - a\ - \b - a\ 

= sign (a — a) (a — a) — sign (b — a) (b — a) 

{sign (b — a) (a — b) if a ^ (a, 6), 
2cr - (b + a) if er e (a, 6). 

Dividing by (6 — a) concludes the proof. □ 

Lemma 4.10. Let E Ax and tp be defined in Lemma \4--6\ For each (x,t) G fly. 

liminf / E Ax (pdyds 



> [ D_[D + sign(A(uA X )-A(u)) 



-(A(uax) + A(SaxUax)) - A(u) 



| ip dyds 



+ liminf i f (( e (u A x,r Ax ,u) + ( £ (ua x ,t£ x ,u)) (A(u) y f(pdyds, 

where 

C(a, b, c) := signed) - A{c)) - sign' e {A{b) - A{c)), Va, b, c e R. 

Proof. We split the proof into two claims. 
Claim 1. 

E) > C £ (z,T+,u)d z (A(z)-Mu 3+1 )) 2 dz 
1 { U] 

+ 2(Axp J C{z,Tr,u)d z {A{z)-A{u j ^ 1 )fdz 
2 



+ \ (C £ («i , tJ , u ) + C £ («,- , r/ , u)) (A(u) y ) 



Proof of Claim 1. Let 



T~ := T^fj ' ^{z,u){A{z)-A{u j - 1 ))dz 

1 [ u i 1 
~ Ax ip"(z,u)dzA(u) y + -dyip' E (u,u j )A(u) y . 
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We start by rewriting the first term as follows: 
1 f Uj 

^(z,u)(A(z)-A(u j . 1 ))dz 

sign' £ (A(z) - A(u))d z (A(z) - ^(%-i)) 2 dz 
signer/) - A{u))d z {A{z) - A( Uj ^)) 2 dz 
( e (z,t- ,u)d z (A(z) - Aiuj^)) 2 dz 



I 



2(A^) 2 Ju^ 
1 



2(Ax) 2 J Uj _^ 
1 



2(Ax)2 J Uj i 
l -s^ e {A{ri)-A{u)){D^A{u 3 )f 



( £ (z,t- ,u)d z (A(z) - A( % -_x)) 2 dz. 



Concerning the second term in the definition of T~ , Lemma 14.81 gives 
ip"(z,u)dzA(u) y = -D_sign e (A(uj) - A(u))A(u) y 

i 

= -signer") - A(u))D_A( Uj )A(u) y . 
For the last term we simply add and subtract to obtain 



f 

Ax 



^d y tl)' e {u, Uj )A{u) y = isign^A^. ) - A(u))(A(u) y ) 2 

^C E {u ,T-,u){A{u) y ) 2 . 



Hence 



T- = Ugn' e (A(r-) - A(u))(D_A( Uj ) - A{u) y ) 2 
1 



2(Ax) 2 J u ._ ± 
+ \( e (u J ,Tr, u )(A(u) v ) 2 . 



C, £ (z,T7u)d z {A{z) - A{uj.{)Y dz 



Define 



1 r Ui 



(Ax) 



in 4 1 



— / ^(z,u)dzA(u)y + -dyip' e (u,u :j )A(u)y. 



AxJuj+r ^ ■ . — 2 
Using the same strategy as above we arrive at 

T+ = isigni(A(r+) - A(u)) (D + A(u 3 ) - A(u) y f 

3+1 



2(Ax) 2 y 



+ ic e (%,r+, W )(A( M ),) 2 . 

Note that E £ = T~ + T + , so Claim 1 follows by removing the non- negative 
on the right hand side. 
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Claim 2. Suppose that x € Ij. Then 



lim inf 



eio 2(Ax) 2 J nT 



C(z,t 3 ,u) — (A{z)-A{u^ 1 )) 2 dz 



+ / C^r+,u)-(A(z)~A(u H1 )) 2 dz 

Ju j + l az 



D- L»+sign(A(u i ) - A(u)) 



Proof of Claim 2. Let 
1 



{A( Uj ) + A(u j+1 )) - A{u) 



ip dyds 

| <p dyds. 
(4.4) 



Ti- 
Kl := 



1 



2(A^) 2 J Uj _, 
1 



2(A^) 2 J Uj+1 
1 



S ign' £ (A(z) - A(u))-(A(z) - A{u^)) 2 dz, 
sign' s (A(rr) - A(u))-^(A(z) - A(u,^)) 2 dz, 
sign^(z) - A(u))^-(A(z) - A{u J+1 )) 2 dz, 
signer/) - A(u))4-(A(z) - A(u J+1 )) 2 dz, 



2(Ax) 2 J Uj+1 " ^ ^ ' dz' 
and note that the left-hand side of (|4.4p may be written 



lim inf / ((Tl - K[) + (Tf - iff)) p dyds. 



Let us rewrite Tf as follows: 
1 



Tf 



(Ax) 2 
i 



(Ax)^ 



i 



si gn ;(A(z) - i4(«))A'(z)(i4(«) - dz 
sign^(yl(z) - A(u))A'(z)(,4(u) - ^(^_i)) 
sig<(A(z) - A(u))A'{z)(A(z) - 



(A*) 2 



where 



JJ? := rr 

1 Ax 2 



sign e (A(z)-A(u))(A(z)-A(u)) 

rf - A(u)\ e dz 



1 dz 

Concerning K\ , we apply Lemma 14.81 to obtain 
1 



K\ 



2(Ax)2 J Uj _ t 



sig<(A(r-) - A(u))d z (A(z) - Aiu,^)) 2 dz 



1 



= -signer) - A( W ))( J D„A(u,)) 2 d* 



D_ S ign e (A( % )-A(n)) J D_A(^)- 
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It now follows that 

Tf - Kf = --^D.sign^Aiuj) - A(uj) 



(A( Uj ) + A(u^)) - A(u) 



R\. 



Performing the same type of computations as above yields 



where 



Kl = ^£»+sign £ (A( % ) - A(u)) 



(A(u j+1 ) + A( Uj )) - A(u) 



2- 



m ■■= 



Ax 2 



f Uj+1 d 
^ -\A(z)-A(u)\ s dz 



S ign e (A(z) - A(u))(A(z) - A( U )) 



z=u j + 1 



Next, observe that 



R{ = sign e (A(z) - A{u)){A{z) A{u)) \A(z) - A(u)\ t 

Ax 

-Z — Uj-i 

so lim 6 ^o i?f = 0. The same considerations apply to i?| so lim e ^o R% — also. 
Claim 2 follows from an application of the dominated convergence theorem. Finally, 
combining Claim 1 and Claim 2 finishes the proof. □ 

4.1. Estimates. The purpose of this section is to find bounds on the "unwanted" 
terms in inequality (|4.2p and Lemma f4. 101 Throughout this section the notation is 
the one given in Lemma FLTjI We let C denote a generic constant. By constant it is 
meant that it does not depend on the "small" variables but it might depend on T 
and the initial data. For any set A, let lyt denote its characteristic function. 
For future reference we collect some standard estimates in a lemma. 



Lemma 4.11. Let ip be defined in Lemma \4-6\ Th 

Qk 

dx 



T^<p{x,t,y,s) 



en 



- Zk+i 1 {\x~y\<r}{x,y)p ro (t - s). 



Recall that S a ip{x, t, y, s) = ip(x + cr, t, y, s). If \a\ < Aa; then 



Qk 



\\p [k) \\r~ 

- W) ^fc+l l {\x-y\<r+Ax}{x,y)Pr {t ~ «)■ 



Considering the difference quotient applied to to r we have 



\D + co r (x-y)\<^-f^l 



{\x—y\<r-\-Ax} 



Proof. Note that 



1 / x \ 

QZ- k ^{x) = -^fl™ (-) 



Since supp(p) C [—1,1] we have 

Qk 



dx k 



u r (x) 



^ r k+l M\x\<r}{?), 



which proves the first statement. 

Consider the second statement. If \x — y\ > r + Ax, then 

\x + a — y\ > \x — y\ — \<j\ > r + Ax — Aa; = r, 

so it follows that l{| I+ff _ 9 |< r }(3;,i/) < ^-{\ x -y\<r+Ax}( x i y)\ this proves the second 
statement. 
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To prove the last statement, recall that 

uj r {x + Ax) — cu r {x) 
D + uj r (x) = — . 

If |x| > r + Ax then uj r (x + Ax) = uj r (x) = 0, so supp(D + (ui r )) C [—r — Ax, r + Ax]. 
By the mean value theorem and the fact that j|w',|| L=0 = ||p'j| L ^ r~ 2 wc get 

\uj r (x + Ax) — uj r (x)\ < — ^2°° 
The last statement follows from this. □ 
Estimate 4.1. 



sign (uax - u) {f{u Ax ) - f(u)) {D + ip + ip y ) dX 
Proof. Let 

fi := j j ^ sign (u Ax - u) {f{u Ax ) - f(u)) {D+p + cp y ) dX. 
First note that 

D+tp + ip y — D + ip - ip x . 

We claim that 

1 f Ax 

(D+ip - ip x ) (x, t, y, s) = — J {Ax - a)(p xx {x + a, t, y, s) da. (4.5) 

Hence 

P = ~r— [ [ [ sign e (A(uax) - A(u)) (f(u Ax ) - f{u)) {Ax - a)S a tp xx da dX. 
Ax JJ U 2 t J 

We can write 

sign {u Ax - u) {f(u Ax ) - /(it)) (x, t, y, s) 

= Y] sign (uj - u) {f{u 3 ) -/(«)) (t, y, s) {x). 

Using summation by parts 
1 



Aa; 



As- 
sign {u Ax - u) {f{u Ax ) - f(u)) {Ax - a)S a ip xx dadx 

1 f Ax f 
= — / V"@j / t {Ij} {x){Ax - a)S a ip xx dxda 

Ax JO j JM 

= — — / y^6j / (f xx {x + a,t,y, s)dx{Ax-a)da 
Ax Jo j Jlj 

/■Ai 

= / Tl^i ( D -S*<Px,j+i/2) {Ax -a) da 
Jo j 

/•Ax 

= - ^ -P+% / S a (p X j +1 / 2 (Ax - a) da, 
j J ° 

where S a <p x ,j+i/2 = <Px(xj+i/2 + &,t,y, s). By Lemma [4.111 we have 
\ip x {x + a,t,y,s)\ < c\t { i x _ y \< r+Ax} (x,y)p ra (t - s). 
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Hence 



f Ax Ax 2 

J o S ff (p Xtj+1 / a (Ax - a) da <C—t^ Xj+i/2 _ yl < r+Ax y(y)p ro (t-s). 



Now 

Therefore 



^ + e J |<||/|| Lip |c 



A.r 



— / / sign (u Ax - u) (f(uA X ) ~ /(«)) (Ax - a)ip xx dadx 



rAx 

<y2\D+Qj\\ S a <p xJ+1/2 (Ax - (7 



)da 



Ax 



< C H/ILpE \ D + U '\ — \\ Xj+1/2 -y\<r + A x }{y)Pro{t *)■ 



It follows by the above and Lemma 157X1 that 



< CAx 1 



Ax r T 



Jo 



= G- 



Ax 



\uax{x + Ax,t) — UAx{x,i)\ dxdt 

n T 

This concludes the proof. 
Estimate 4.2. 

\A(u Ax ) - A(u)\ [D-D+v + (D+ + D-)<p y + y yy ) dX 



□ 



< C- 



.Ax A 



1 



Ax 



Proof. Since ip xx + 2<p xy + tp yy — it follows that 

D-D+ip + (D+ + D_)<p v + ip yy = {D _ D +ip - ip xx ) + ( (D + + D_)ip - 2(p x ) y 
Thus 



\A(uax) - A(u)\ {D-D+ip + (D + + D_)ip y + (p yy ) dX 



\A{uax) - A(u)\ {D^D + ip - <p xx ) dX 



+ J J \A{uax) - A{u)\ ((£>+ + DJ)<p - 2<p x ) y dX 

= ■ Cl + C2- 

Consider the term (i. We use the same strategy as in Estimate 14.11 Writing 
/i(cr) = ip(x + a, i, y, s), a Taylor expansion gives 

H(z) - M (0) = V(0) + ^V'(O) + ^V 3) (0) ~\[^~ ^)V (4) (^) da. 

Using this, we get 

//(Ax) - 2/i(0) + fj,(-Ax) - Ax 2 n"{0) 
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i r Ax i r° 

= — / (a- Ax) 3 /j, {4) (a) da + - (a + Ax)V 4) (a) da. 
6 Jo 6 J -Ax 

It follows that 

1 f Ax d 4 
D + D-<p - <p xx = - — 2 / (a- Ax) 3 — — <p(x + cr,t, y, s)da 
6Ax Jo ox A 

+ eh? Lj a + ^h v{x + a > *• 2/1 s) dG - 

Splitting Ci according to this equality we get 



Ci = / 1 2 |A(uAx) ~ A(u)| (£-.0+^ - Vxx ) dX 

1 f f f Ax d 4 

\A(u Ax ) - A(u)\ (cr - Ax) 3 — -jip(x + a, t, y, s) da dX 



6Ax 2 JJn^Jo dx 4 

1 f f f° d 4 

\A(u Ax ) - A(u)\ (a + Ax) 3 —<p(x + a,t,y,s) da dX 



6Ax 2 JJ n 2 T J-Ax'"' ~~' dx 4 

=■ Cl,l +Cl,2- 

We also have that 

\A(u Ax ) - A(u)\ (x, t, y, s) = ^2 \A{ Uj ) - A(u)\ (t,y,s) l {Ij} (x). 

3 ' % ' 

Now consider £1.1, 

/ \A(u Ax ) ~ A(u)\ (a - Ax) 3 — S a tpdxda 
Jr ox 

/•Ax r Qi 

= -22\A( Uj )-A(u)\{t,y,s) ^ (a - Ax) 3 J^t {I]} (x) — S a (pdxda 

/■Ai 

= -Ax / (a - Ax) 3 ®jD_<p° XX!j+1/2 da 

/•Ax 

= Ax D+$j / (a - Ax) 3 S a (p X xx,j+x/2da, 
Jo 



where 



Q3 

Safxxx,j+l/2{t, V, s) = -g^(p(x j+ i/2 + a, t, y, s). 



Now we use Lemma [4.111 to estimate this term as follows: 
1 f f f Ax d 4 

iCi,ii= te^JJ J ^ UA ^~ A ^^ a ~ Ax ^d^^ x + ^^y^^ d(7dx 

1 f f T f Ax 
6Ax J J ^2 D+ ® j J - Ax) 3 S a iPxxx,j+i/2da dtdyds 



r + At r / r Ax \ 

<C , , / \D+A(u Ax )\( (a ~ Ax) 3 da) dxdt 
Ax r 4 Ju T Vo ' 



Ax'r 4 Jn T 

^ 2 r_±Ax =c A£ ( i + Ax 

r 



where we have used that \A(iiAx{-,t))\ BV ^ is bounded independently of Ax 7 t,rj 
by Lemma |3. II The term is estimated in a similar way. 
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Now consider £2- Again, let /x(cr) = ip(x + a, t, y, s). Then 

(D+ + D_)tp -2<p x = ^- Ifi(Ax) - fx(-Ax) - 2Ax M '(0)] . 
lax 

By a Taylor expansion 

H(z) M (0) = V(0) + ^V(O) + 5 jfV - z) 2 ^(a) d<r. 
Puting 2 = ±Ax and subtracting the corresponding equations we obtain 
{D + + D_ )(p - 2(p x = J (a- Ax) 2 -^ip(x + a, t, y, s) da 



+ 



2Aa 



Q3 

(a + Ax) 2 -g^(p(x + a, i, y, s) da. 



Ax 



We may split (2 into the two terms 
C2 



1 

2 Ax JJ n 2 J 



d 3 d 



\A(uAx) - A(u)\ ( a - & x ) 2 g x 3 g y ^ X + CT '^y' S ) da dX 



1 ff f° d 3 d 

^ J L l-A x lA(UAx) " A{U) 1 {a + ft? ^ + *' *' " S) ^ 



2A 

= : C2,l + C2,2- 

Performing integration by parts, £2.1 becomes 



1 r r r Ax gs 

^~xJJ^J s ' 1 S n (A(uAx)~A(u))A(u) y (a-Ax) 2 —Lp(x + a,t,y,s)dadX 



2A 

Thus, by Lemma |4.11[ 
lCa.il < 



2A.t 

r + Ax 



\A(U) V 



< TC 



r 4 Ax 



f Ax d 3 
J (a- Ax) 2 T^ip(x + a,t, y, s) da 

[ (a - Ax) 2 da j f \A(u) y \dyds 
Jo J Jn T 



dX 



< C 



Ax 2 



A., 



as \A(u(-,s))\ BV(R) < \A(u°(-))\ BV(m for all s. The same estimate holds for ( 2 ,2- 



□ 



Estimate 4.3. 



sign [z — u) F 2 (z) dz D + ip dX 



.Ax 



< C — 1 + — 



Ax 



Proof. By definition F 2 is bounded. Hence, 



/ sign (z-u)F^(z)dz <\\F 2 \\ Up Ax\D +Uj \. 

J Uj 

Note that \ua x (' ji^gy,-^ is bounded independently of Ax,t,r] by Lemma l3.lt so 



we may apply Lemma 14.111 to obtain the result. 
Next, we consider the terms from Lemma [4.101 



□ 
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Estimate 4.4. 

D- (D+sign^^) - A(u)) 



{A( Uj ) + A(u 3+1 )) - A(u) 



ipdX 



A.! 



> -C(l + r + Ax) — 1 + 



Aa 



Proof. Let us first show that 
D + aiga(A(uj) - A(u)) 



First note that 



(A( Uj ) + A(uj +1 )) - A(u) 



< AxL» + sign(A(w i ) - A{u))D + (A( Uj )). (4.6) 



D+siga(A(uj) - A(u)) 
2 

= ^sign(A( Uj -) - A(u J+1 ))l {A{u]eint{AiUjhA{Uj+l)} , 

so the left-hand side of (|4.6[) is zero whenever ^ int(A(itj), j4(uj+i)). Second, 
if c £ int(a, b), then it follows that 



(a + b) - c 



1 



(|6-c| + |o-c|) < |fc-o|. 



Since z i-> sign (A(z) — A(u)) is increasing, the right-hand side is positive. This 
proves (|4.6I) . 

Performing integration by parts we obtain 

1 



' D- ( D + sign(A(uj) - A{u)) 
D + sign(A( Uj ) - A(u)) 



-(A( Uj ) + A(u j+1 )) - A(u) 



-{A( Uj ) + A(u j+1 )) - A(u) 



< Ax I J D + sign(A(uj) - A(u))D + (A(uj)) \D + (p\ dX. 



<pdX 
D + ip\ dX 



Using integration by parts for difference quotients and the Leibniz rule for difference 
quotients, we obtain 

J I ^ D + sign e {A{uAx) - A(u))D + A(u Ax ) \D + <p\ dX 

= - JJ^ sign £ (^(u Aa: ) - A{u))D + A(uax)D- \D+tp\ dX 

sign E (A(u Ax ) - A(u))D-D + A(u Ax ) \D_ip\ dX 



= :Ci+C2- 

To estimate £i we first observe that D_ < \D+D-<p\. Furthermore, when 

proving Estimate 14.21 we established that 



D + D_tp(x,t,y,s) 
= ip xx (x,t,y,s) 



6Ax 



Ax Qi 

(a - Axf-^-j(p(x + a, t, y, s) da 



Lax {(J + * xf h v{x + a > y > s) da - 
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By Lemma ETTTl 



f±Ai 



d 4 

(cr =p Axf-^ip{x + a, t, y, s) da 



(Ax) 4 

< C—^—^{\x-y\<r+Ax}(x,y)p ro (t - s) . 

Using Lemma 14. 1 II once more, the above implies that 



<-j 


\<Pxx 


dyds - 




<A 


1 + 


Ax 2 \ 






r 4 J 



Ax 



Therefore, 



ICi 



< 



JJ^ sign e {A(uAx) - A{u))D + A{u Ax )D+ \D-tp\ dX 

I \D + D-<p\ dyds) dxdt 
Jn T ' 

1 + — ) / \D + A{u Ax )\ dxdt. 
r J Jn T 



D+A(u Ax )\ 

2 



Recall that \A(u Ax (-, t))\ BV ^ is bounded independently of Ax, t, r\ by Lemma |3~T1 
Concerning £2 we have 

IC2I = 



< 



sign £ (A(u A;c ) - A(u)) (D^D+A(u Ax )) \D-<p\ dX 

J J \D_D+A(u Ax )\ \D-tp\ dX 
\x f 

— / \D-D+A(u Ax )\ dxdt. 



< C— 

r 

Note that it follows from (|3.1|) and Lemma [3.21 that \\D^D + A(u Ax (-,t))\\ Ll ^ is 
bounded independently of Ax, t, r\. Hence, 

Ax J J D+sign £ (A(u A x) - A(u))D + A(u Ax ) \D+tp\ dX 
<L(|Ci| + |C 2 |) 



Ax 



< C(l + r + Ax)^- 



Ax 



□ 



Estimate 4.5. 



Um inf 



(C e {uAx,T Ax ,u) +C(u Ax ,t a ~ x ,u)) (A(u) y ) 2 ipdyds\ dxdt 



eio 7 nT 2 

„ ( Ax Ax Ax 
>-C [ — + — + — 

Proof. Set 

TZ] := (C £ K,r-,u) + C( Uj ,t+ 7 u)) (A(u) y ) 2 

= (sign^( % ) - A{u)) - signer-) - A{u))){A{u) v f 
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+ ( S ign' £ (A( % ) - A(u)) - sig<(A(r+) - A{u))){A{u) y f , 

and lZ e Ax (x, t, y, s) = lZj(y, t, s) for x S Ij. Note that the term we want to estimate 
may be written 

liminf / TZ% x {x,t,y, s)tp(x,t,y, s) dyds 

= X] 1 ™ 1 ^ yj ^(t^, s)Lp{x,t,y, s) dyds^j l {Ij} (x). 
Let us define an entropy function by 

:= sign £ (A(07) - A(u)) - 2sign e (A(^) - A(u)) + sign £ (A(0+) - A(u)). 

Recall that ffj = 6 j(u) , so the above function is not as explicit as it appears. 
However, by Lemma l4.9l we are able to obtain an explicit expression for the limit as 
e — > 0. To simplify the notation wc write ^(u) for d u ^ s {u, Uj—i, Uj, Let 
us also define the entropy flux functions 

That is (*&e,ji 3c,j! ®e,j) is an entropy-entropy flux triple. 
Multliplying equation (jTTTJ) by ^f' e • (it) yields 

+ S eii (u)„ = $ e j(u)„„ - d v % tj {u)A(u) y . 
By Lemma T4. 81 wc sec that 

dy%j(u)A(u) y = [sign e (A(^) - A{u)) y - 2sign £ (A(u j ) - A(«)) w 

+ sign e (A(6f)-A(u)) y ]A(u) y 

= - signer) - A(u)) - 2sig4(A( % ) - A(«)) 



+ S ign' s (A(rf)-A(u)) {A{u) y f 



It follows that we can write 

TZjtp dyds = / $F Stj {u)ip s + E Stj (u)ip y + (u)ip yy dyds 
=:T[+Ti+T£. 

Let us consider the three terms separately 
By Lemma l4~9l 



(4.7) 



(4.8) 



lim ^ j (u) = lim / ■ (z) dz = / lim ■ (z) dz 

e;0 ' eiO J Uj ' J J u . eiO 

S 

sign(A(z) - A(uj)) - s gj _! (A(z)) dz 
+ / sign(A(z) - A(« 3 -)) - sgj (A(z)) dz, 



where sg^ (cr) := sg (A(Mj . ) A(Uj+l)) (cr). Again by Lemma the mapping 
z ^ sign(A(z) - A(uj)) - sgj-j (A(z)) 
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has support in mt(uj, Uj—i). Similar considerations apply to the second term. 
Hence 



lim \& e 

ej.0 



< 



sign(A(z) - A(uj)) - sg, ! (A(z)) d« 



3ign(A(z) - Afuj)) - sgj (A(z)) rfz 



< 2 |uj — + 2 |ttj+i - Uj| 

By the same type of reasoning we obtain the bound 



lim Sy(«) 



< 



3 ign(A(z) - A(«,)) - sg,.! (A(z))] f'(z) dz 



u j+1 



[sign(A(z) - A( Uj )) - s gj (A(z))] f(z) dz 

< 2||/'|| ioo (|tij + \u j+ i -Uj\). 

Concerning we use substitution and the explicit expression given in Lemma 
This leads to 



lim<£>^(u) 



(-sg,-,! (A(z)) + 2sign(,4(z) - A(uj)) - sgj (A(z))) A\z)dz 



A(u) 



< 



A( Uj ) 
A(«j_i) 



sgj-i (o-) + 2 sign(er - A(tij)) - sgj (cr) dcr 
/ sign(cr - A(/uj)) - sg, J (cr) da 

J A(u 3 ) 

r-A(u j + 1 ) 

/ sign(cr - A{uj)) - sg - (ct) da 

JA(ui) 



'A(uj) 

< \A{uj) - A(uj-x)\ + \A{u j+1 ) - A( Uj )\ 



Let us return to equation (|4.8[) . By the dominated convergence theorem and the 
above computations 



\y s \dyds<C{\D„u 3 \ + \D + u 3 \) 



limTf 

ej.0 


< 


lim ^ 

£^0 J 


limTl 

e|0 


< 


lim 


limX? 

ej.0 J 


< 


lim 

£^0 J 



L°° Jrir 



.Arc 



|^|<%ds<<7— (iD.Uj-l + ID+ttj-l), 
Ax 

\<p yy \ dyds < C— (\D_A( Uj )\ + \D + A{u J )\) 



Hence 



/ V^liminf I / TZ £ {t 1 y,s)ip{x 1 t,y 1 s)dyds\ t^j j y(x)dxdt 

J YIt j e ^ \J IIt / 

{Ax Ax\ f 
> -C + / 1-D-UAzl + |£>+uaz| ctecft 

Ax f 

- C— / |D_A(w Ax )| + |D + i4(«Ax)| dxdt. 
r Jn T 

The desired estimate now follows from the uniform bounds in Lemma 13. 11 □ 
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4.2. Proof of Theorem 14. 11 Let us now combine the previous results to conclude 
the proof of Theorem 14.11 We begin by stating a rather standard lemma. 

Lemma 4.12. Set 

n{t) := 



\u A x(x,t) - u(y,s)\oj r (x - y)p ro (t - s)dydsdx. 
Let t > ro, and denote by L c the Lipschitz constant of t ^ \[u(-,t)\\ L \(^y Then 

n(t) - \\uAx(;t) - u{-,t)\\ LHm < \u{-,t)\ BV(R) r + L c r . 
Proof. By the reverse triangle inequality, 
«(*) - \\u Ax (-,t) -u(-,t)\\ L1Q 



< 



< 



\u(y, s) — u(x, t)\ uj r (x — y)p ro (t — s) dydsdx 



r Jn T 

T 



\u(y, s) - u(t, y)\dy p r<) (t - s) ds 



\u(t, y) — u(x, t)\tu r (x — y)dydx 
< L c r + \u(-,t)\ BV(R) r. 



□ 



Proof of Theorem \4.1\ Our starting point is Lemma 14.61 Let A{a) = A(a) + na, 
where A is the original degenerate diffusion function. Let 



sign (u Ax - u) (f(u Ax ) ~ f(u)) (D + ip + <p y ) dX 



sign (z — u) F^z) dz J D + ipdX 



+ JJ r> \A u Ax) - A(u)\ (D-D+ip + (D + + D_)ip y + ip yy ) dX. 
By Estimate 14.11 Estimate 14.21 and Estimate 14.31 it follows that 



131 <C 



Ax 



1 



Ax 



1 



Ax 



--:E X . 



(4.9) 



Furthermore, by Lemma 14.71 Lemma 14.101 Estimate 14.41 and Estimate 14.51 it 
follows that 

liminf// E Ax ipdX > -C{l + r + Ax)—^- ( 1 H ) - C =: -E 2 . (4.10) 

eiO J Jxi 2 T r V f J r 

Applying the estimates f)4.9[) and (|4.10[) . the inequality (|4.2[) becomes 
\uax -u\p a (t- t)u t (x - y)p r „(t - s)dX 



< 



\uax - u\ p a (t - v)uj t (x - y)p ro {t - s) dX + Ei + E 2 . 



Note that both E\ and E 2 are independent of a. Thus, we can send a to zero, 
arriving at 

k{t) < k{v) + E 1 +E 2l 
where k is defined as in Lemma 14.121 By Lemma 14.121 it follows that 
\\uax(-,t)~u(-,t)\\ lhr) 
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< \\uax(;v) -u(-,u)\\ L i (M) +2 {l c r + |u°| w(R) rj + E 1 + E 2 . 

Recall that we had to pick v > ro. Denote by Ld the L 1 Lipschitz constant of 
t 1 — ^ t). By the triangle inequality 

\\UAx(-,v) -U(; f)|| L i (K) 

< \\uax(; v) - wL|Li (K) + - "°Li (8) + IK - «(•, ^)ILi (R) 



This means that 



\\uax(-,t) - u{-, t)|| l1(r) < - u°|| L1(R) + (L c + L d ) v 

+ 2(L c r +|u°| BV(a) r)+£ 1+ £ 2 . 

Choose r 3 = Tq = Aa; and ^ = 2ro. Then there exists a constant C such that 

||«a*(-,t) -u(;t)\\ l1Qsl) < +CAx*. 

Now recall that A(a) = A{a) + rj<r and so we need to send r\ to zero to finish the 
proof. If u n is the classical solution of the regularized equation and u is the entropy 
solution of the non-regularized equation, then it is well known that t) — > u(-,t) 
in L 1 (M) as rj — > (see section^. Concerning the scheme one may prove continuous 
dependence in I 1 on 77 using Gronwall's inequality. Hence, we can also send n to 
zero in the scheme. This finishes the proof of Theorem 14.11 □ 



5. Implicit difference schemes 

In this section we show that the arguments presented in the previous sections 
carry through for implicit schemes. Fix a time step At > 0. We consider implicit 
difference schemes of the form 

D l _u n 3 + D-F(v%, u? +1 ) = D-D+ A(v$) n > 1, 3 G Z, (5.1) 

where 

D < = 3 3 — . 
" 3 At 

Let t n = nAt and Xj = jAx. We define the grid cells 

7™ = [xj-x/2,Xj +1 / 2 ) x (t n -i,t n ] for n > and j G Z. 

The piecewise constant approximation is defined for all (x,t) <G R x (— Ai,T] by 

UA(ar,t) = u" for (x,t) G Ij\ (5.2) 

The domain is chosen so that D^ua is defined for all (x,t) E R x (0,T). For the 
existence of a unique solution u" to the nonlinear equation (|5.1[) and the convergence 
of ua to an entropy solution, see |llj . 
We now state the main theorem. 

Theorem 5.1. Let u be the entropy solution to (|l.ip . and let ua be defined via 
u™ by (|5.2p . where u™ solves (|5.ip . // u° satisfies the same assumptions as in 



Theorem \4-l[ then for all sufficiently small Ax and At, and for all n G N such that 

tn G [0,T], 

\\u A (;t n ) - u(-,t n )\\ L1{R} < \\u° A - u°\\ L1(R) + C (Ax 1/3 + Ai 1 / 2 ) , 
where the constant Ct depends on uq, A, /, T , but not on Ax, At. 

To prove this theorem we will follow step-by-step the proof of Theorem 14.11 and 
present the details whenever there is a significant difference between the two cases. 
Thanks to [TTJ Lemma 2.4], we have the following L 1 Lipschitz continuity result: 
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Lemma 5.1. Let m and n be two non-negative integers. Then 

|| u a('i t n ) — Ua('j ^itOIIl^R) — \t n — t m \ , 
where L d = \F(u°,u° +1 ) - D+A(u*)\ BV . 

Next, let us prove an implicit version of Lemma 14.31 
Lemma 5.2. Let u™ be the solution to (|5.1[) ■ Then for all c € R, 

D*_Mv?>c) + D_Q c (u],u] +1 ) - D_D + \A(v%) - A(c)\ 

<-7a^2 r ^(z,c)(A(z)-A(u"))dz 



3+1 



where Q c (u,v) is defined in Lemma \4-3\ 
Proof. From (|5.1[) it follows that 

tt(tf,c)Dtu] + rl>Uu?,c)D-F(u?,vZ +1 ) = ^ e {v7],c)D_D + A{tf). 
Apply Lemma 14.21 with g(cr) = a, a = uj, and b = u"^ 1 to obtain 

^(u7, C )D*_u7 = D*_VsK",c) - ^ / ^ ^'(*,c)(* - u?" 1 ) 

The remaining part of the proof follows exactly as in the proof of Lemma 14.31 □ 
Let us define the time shift operator 

S t M a(t)=a(t + M), 

for any function a = o~(t). 

Lemma 5.3. Suppose A' > 0. Let ua = u/±(x,t) be defined by (|5.2p . and let 
u = u(y, s) be the classical solution of (jTTTJ) . Let ip(t) := ln^ T \\(i) and define 

tp(x,t,y,s) = ip(t)u r (x - y)p ro {t - s), 

where ui ri p ro ,v, r are chosen as in Lemma \4.6\ Then 

| it A - u\ 5 At (t - v)u r p ro dX 

n'l 



\u A - u\ S At ^u; r (D + p ro - d t p ro )dX 

+ At I I \UA — U\ D+ilUJrdsPra dX 



+ J J a sign (u A - u) (f(u A ) - f(u)) (D + tp + tp y ) dX 

r r ( / > ^ AxtiA \ 

+ / / / sign(z - u)F^(z)dz) D+tpdX 

J Ju% \Ju Ax J 

+ JJ 2 \Mua) ~ A(u)\ {D_D +V + (D+ + D-)<p y + (p m ) dX 

> // \ua — u\ 5 At (t — T)ui r p ro dX + ]hnixd // E A ipdX, 
J Jul £ W J Jul 
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where 

<*At(*) = A7 1 {[-At,0)}(*) ! 
andE A (x,t,y,s) = E'[u] ,u? +1 )(y, s) for (x,t) e J?. 
Proof. As in Lemma 14.51 we obtain by Lemma 15.21 the inequality 

- {dl + 0„(I>_ + D+) + D_i? + )|^(u7) - A(«)| e ) < 

where Ej n := u™, is defined in Lemma 14751 Let us multiply by tp 

and integrate over LT^. Integration by parts for difference quotients and ordinary 
integration by parts gives 

1pe(uA,u)D t + if + l/) e (u,UA)tPs dX 

2 
T 

+ JJ ^ q E (u,u A )tp y + Q u (u A , SAxUA)D + ipdX 
+ JJ n2 \A(u A ) - A(u)\ e {<p yy + (£>_ + D + ) Vv + D_D +V ) dX 
> J J E%ipdX. 

Consider the first term on the left. Let e tend to zero as in the proof of Lemma l4~6l 
Using the Leibniz rule for difference quotients and adding and subtracting we obtain 



D\ip = S At ll)UJ r D\p ro + D+V^rPro- 



Furthermore, 
Hence, 



ifi s = -S At lpOJ r d t Pr + AtD+tpUrdsPro- 



u A - u\ (D+tp + ip s ) dX 

\ua - u\ SAt^T-C-D+Pm - dtPr ) dX 

+ At \u A — u\ D+lpbJrdsPrg dX + \llA — u\ D+lpUlrPrg dX . 

J Jul J Jul 

Finally, we use that 

D%i> = 8 At {t-v)-5 At {t-T). (5.3) 
The lemma now follows, as in the proof of Lcmma l4.6l by letting e tend to zero. □ 

Comparing the terms in Lemma 14.61 with the terms in Lemma 15.31 we recognize 
all but two terms. 

Estimate 5.1. 

|U - Ua| S^lpOJriD+PrQ - d t p ra ) dX 

Proof. To show this we use a Taylor expansion: 

p ro (t + At-s)- p ra (t - s) 
r t+At 





At\ 











J -Q^PnX z - s ) dz 
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dz 



ft+At Q2 

p ro (z - s)At - J ~d£2 Pr °( Z ~ ~( t + A< )) dz - 



It follows that 

V 



D + p ro - d t pr = 

Integration by parts yields 



t+At q2 

At J t dz^ 



p ro {z-s)(z-(t + At))dz. (5.4) 



1 

At 



t+At fT 



\u - u A \^——p r (z - s)(z — (t + At)) dsdz 
oz 1 



At Jt 



t+At r T 



d d 

\u - u A \ E — —p r „(z - s)(z- (t + At)) dsdz 



i r t+At r T d 

At J J s1 S u e( u - u &) u s-g^Pr {z - s)(z - (t + At)) dsdz. 



Since 



o Pro\ z s ) 

oz r oz \ r 

and p ro has support in [-ro.ro], it follows that 



19 ( z — s\ 1 , I z — s 



1 

At 



t+At r T 



< 



c 

'rjAtj, 



d 2 

\u(y, s) - u A (x,t)\ e —p ro (z - s)(z - (t + At)) dsdz 

/t+At 
l{|z- s |<r }|2 - (* + At) | dzds 



At f T 

<C— |u S |l{|t-s|<r +At} d,S. 

r o Jo 



Multiply the above inequality by ip(t)ui r (x — y) and integrate in x,y,t. From the 
resulting inequality and (|5.4[) . we arrive at the estimate 



\u - u A \ £ S At ipui r (D + p ro - d t p ro )dX 



< c- 



Y [[ |u s |5^ t ^w r l { | 4 _ s |< ro+At j dX < C— (l + — ) ||M s || L i ( n T ). 
o J Jti*, r o \ r o J 

Since ||u s ( - , s)||l 1 (e) is uniformly bounded on [0,T], the estimate follows from the 



dominated convergence theorem. 
Estimate 5.2. 



< C- 



At \ll — U A \ D+lpUJrdsPrg dX 

J Jill 

Proof. Integration by parts yields 

J J ^ U^D+lpUJrdsPrgdX < J J ^ \u s \\D\ 

Because of (|5.3|) and since 



At 



ip U! r p r() dX. 



\\Pro\\L<*> ^ 



it follows that 



\D\tp\ u) r p ro dX < 2 ^|j L °° || Ug || l i {YIt} 



ro 



□ 
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Proof of Theorem \5.1\ We start out from Lcmma. [5.3l with A(o~) = A{a)+rja^ where 
A is the original degenerate diffusion function. By Estimate 15.11 and Estimate 15.21 



\u — ua\ S At tpuj r (D + p ro - d t pr ) dX 



At// \u-u&\D t ,'il)U} r d s p ro dX < C— ( 1 + — I =: E 3 . (5.5) 
JJu% r a V r o 1 



Since all the estimates from Section 14.11 apply, we obtain 
|ua — u\ S At (t ~ r)ui r (x — y)p ro (t — s) dX 



At 



^JJ 2 \uA X -u\8* t (t-v)uj r (x-y)p ro (t-s)dX + E 1 + E 2 +E 3 , 

where E\ and E 2 are defined respectively in (|4.9[) and (|4. 10[) . 

Let us make the simplifying assumption that v = t m and r = t n for some 
m, n € N. Then the above inequality rewrites as 

K(t n ) < n(t m ) + E l + E 2 + E 3 , 

where 

K (t) = / / \uA{x,t) - u(y, s)\ cu r (x - y)p ra {t - s) dydsdx. 



Applying Lemmas 14.121 and 15.11 and following the reasoning given in the semi- 
discrete case, we arrive at 

\\u A {-,t n ) - u(-,i„)|| L1(R) 

^ || "A « ||il(R) 

+ [L c + L d ) t 

a s2 (, Ax\ 3 Ax Ax „At A At\ 
+ C(l + r + Ax) 2 1 + — — + C — + C — 1 + — 
V r J r r o r o V r o / 



< IK - u °\\mm + C I -ir + + r + r 



'Ax Ax + At 

where is the constant in Lemma 15.11 and L c is the constant from Lemma 14.121 
Minimizing over r and ro, it is straightforward to see that for sufficiently small At, 
the minimum of the last term is dominated by 

C (Ax 1 / 3 + At 1 / 2 ) . 

This proves the theorem. □ 



6. Explicit difference schemes 

In this section we use the techniques developed in the previous section to provide 
a similar result concerning the explicit scheme. Fix a time step At > 0. We consider 
cxplicit schemes of the form 

D f + u] + D_F(u],u^ +1 ) = D_D + A(u]) n > 1, j G Z, (6.1) 



where 
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The relevant a priori estimates and convergence to an entropy solution is proved in 
[12"] under the hypothesis 

l-^(Fl(z)-F^z))-2^A>(iv)>0, V(z, W )eM 2 . (6.2) 
i±x Ax 

Let t n = nAt and xj = jAx. We define the grid cells 

/" = [Xj-i/2,x j+ i/ 2 ) x [t„,t„ + i), for n > and j e Z. 

The piecewise constant approximation is defined for all (x, t) G ITy by 

u A {x,t) = u] for (x,t) € Z™. (6.3) 

Theorem 6.1. Lei it &e the entropy solution to and letu^ be defined by (|6.3[) 

no m™ ; where u™ solves (|6.1[) . Suppose At and Ax are chosen such that (|6.2[) arte? 
£/ie strengthened condition At < CAx s ^ 3 hold. Ifu° satisfies the same assumptions 
as in Theorem \4- 1\ then for all sufficiently small Ax, and for all n £ N such that 

tn € [0,T], 

\\u A (;t n ) - u(;t n )\\ LHK) < \\u% - u°\\ Ll{K) +CAx 1 ' z , 

where the constant Ct depends on A, f,u° , T , but not on Ax. 

We begin by proving the following lemma 
Lemma 6.1. Let be the solution to (|6.ip . Suppose that 

l-^(F{(z)-Fi(z))>0, (6.4) 

for all z S IR. Then for all jgZ and neN, 

D%tP s (u],c)+D_Q%u],u? +1 )-D_D + \A(u?)-A(c)\ £ 
1 

(Ax) 



< -7X32 / C )( A W - A ( U U)) dz 



'3-1 

1 



(Ax) 2 7„ 



'3 + 1 
+ 1 



^,c)(A(z)-A(^ +1 ))<fe 
^'(.z,c)<fej D_D+A{u]), 
where Q c (u,v) is defined by 

ru rv 

Q c (u,v)= / ip' s (z,c)F[(z)dz+ / ip' s (z,c)F2(z)dz, u,v € 



Proof. We divide the proof into two steps. 
Claim 1. Let {m™} be a solution to (|6.ip . Then 

D^u^c) + D_Q c (u^u^ +1 ) - D_D+\A(u]) - A{c)\ 
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where 

1+1 



— / TP>:(z,c)(z-u] +1 )dz 



1 f u i 

3 — 1 

Ail 



(F 2 (z) - F 2 (u] +1 )) - ±(A(z) - A(u] +1 )) 



dz 
dz. 



Proof of Claim 1. By definition (|6.1[) of {it"} it follows that 

V>eK,c) [D*.u7 + D-F{u],u% +1 ) - D-D+A(u])] = 0. 
Let g(z) = z in Lemma |4~21 It follows that 

n+l 

^( ttj " lC )^ = %(^, c ) + - y ^'(z, c )(z-^ +1 )dz 

The remaining terms can be treated as in Lemma 14.31 
Claim 2. Suppose holds. Then 

i+i 



1 



— / ^'(z,c)(z-^ +i )dz 



1 

Ax 
1 

Az; 



3-1 
""+1 



> 



^'(*, C )(F 1 («)-.Fi(uy_ 1 ))d2 ! 

^(z,c)(F 2 (z)-F 2 (u] +1 ))dz 

-1 

^{z 1 c)D^D+A{u r l)dz. 



(6.5) 



Proof of Claim 2. Consider the first term on the left-hand side of (|6.5[) . By 
definition, = u™ - AtD-F(u 7 f, + AtD-D+A(u 7 f), and so 



1 

a7 



^'(z,c)(z-< +1 )^ 



^ / ip"(z,c)(z-u r })dz 



^{z,c)D^F(u , l,u , l +1 )dz 



i/%(z, c)D_£>+A(u") =: Ti + T 2 + T 3 



Note that Ti is positive. Let us split T 2 according to 

D_F(u?,u? +1 ) = £ (Ft(v?) - F 1 {^_ 1 )) + £ (F a («y +1 ) - F a («y)) 



and thus 



To 



Ax 



tf(z,c) (F 1 (u^-F 1 (u^_ 1 )) dz 



+ — I ^(z, c) (F 2 (uj- +1 ) - F 2 (u])) dz. 
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Now, let us split the two other terms appearing in equation (|6.5p : 
5i:=— / ^{z,c)(F 1 {z)-F 1 {u 1 y_ 1 ))dz 



Ax 



3-1 



3-1 

„," + ! 

1 



Ax 



^(z,c)(F 1 (z)-F 1 (u 1 l_ 1 ))dz 



and 



1 /-"i+i 

S 2 :=— I ^ ^(z,c)(F 2 (z) - F 2 («y +1 )) rfz 



^'(^c)(F 2 (z)-F 2 ( U ;Vi))rfz 



"3+1 
1 / >t '" +1 

+ — / 3 C(^c)(F 2 ( Z )-^ 2 « +1 ))^. 



7/ 



Combining the above expressions we obtain 

i 

T 2 +S 1 + S 2 = -— / ^ ^{z,c){F x {z) - Fi(v$)) dz 



1 

Ax 

1 

Aa; 

1 

Ax 



3-1 



^(z,c)(F 2 (z)-F 2 (u]))dz 
^(z^^F^-F^u^dz 
^(z,c)(F 2 (z)-F 2 (u? +1 ))dz. 



j + 1 

The two last terms on the right-hand side are positive as F is monotone. Let 

H(z) = z-^(F l (z)-F 2 (z)). 
Then, by assumption (|6.4I) . 

i r u i 

T 1 +T 2 + S 1 +S 2 >— J ^ ^(z, c) [H(z) - H(u])] dz > 0. 

Adding T3 to both sides proves Claim 2. 
By Claim 2 



1 

> 



Ax z J u 



^(z,c)(A(z)-A(u]_ 1 ))dz 



/n <'(*,c)(A(*)-A(^ +1 ))^ 
?+i 



Aa; 

.+1 

ip^(z,c)D-D + A(u^)dz. 
Combining this with Claim 1 proves the lemma. □ 



ERROR ESTIMATES FOR CONVECTION-DIFFUSION EQUATIONS 



37 



Lemma 6.2. Suppose A' > 0, and (|6.4p applies. Let u A = uj\(x, t) be defined by 
(|6.3[) . and let u = u(y,s) be the classical solution of (jl.lj) . Set ip(t) := 1{[iat)}(^) 
and define 

(p(x,t,y,s) = tp(t)u r (x - y)p ro (t - s), 
where u> r , p ra ,v,T are chosen as in Lemma \4-6] Then 

\ua - u\ 5^ t (t - u)u) r p ro dX 

| MA - u\ S t _ At 4>0J r {D i _pr - dtp r „)dX 

I j ^ |«A _ u l D^lpUJrdsPrg dX 

sign (ua - u) (/(""a) - f(u)) (D + ip + <p y ) dX 



// / sign (z-u)i^(z)dz } D + ipdX 

JJn% \Ju Ax J 



+ / / \A{u A ) - A{u)\ {D_D+ip + (D+ + D^)<p y + Vyy ) dX 



> 1 1 \uA-u\5^ t (t-T)u) r p ro dX 

liminf // E%wdX 
<40 JJiP T 

At L>^.sign 6 (A(w A ) - A(u))D-D+A(uA)<pdX, 



whe 



5 At(*) = -tt1{[o,a*)}(*)> 



At 

andE A (x,t,y,s) = E^[u]{u^_ x , u],u] +1 )(y,s) for (x,t) e I?. 

Proof. By Lemma 16. 11 wc obtain as in Lemma 14.51 the following inequality: 

D^Mu^u) + d s Mu,u]) + d v q £ {u,u n 3 ) + £_q>;>™ +1 ) 
- (d 2 y + d v (D_ + D+) + D-D+)\A(v$) - A(u)\ £ ) 

< -E]. n + (yJ2 tf(z,u) dz^j D_D + A(u]), 

where Ej n :— E £ [u](u™_ 1 ,u™ is defined in Lemma |4"31 Note that 

^(z,u)dz = At D\.siga e (A(u]) - A(u)). 
Integration by parts for difference quotients and ordinary integration by parts gives 
iK, u)Dt_ip + ip e {u, "aVs dX 

+ JJ ^ q e (u,u A )fy + Q u (u A ,SAxUA)D + ipdX 

+ 11 \A{u A ) - A(u)\ £ {<p y y + (£>_ + D + ) Vy + D_D +V ) dX 
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> 



// E A ipdX-At Disign £ (A(it At ) - A(u))D_D + A(u A )<pdX. 
J Jul J Jul 



Consider the first term on the left-hand side. Let e tend to zero as in the proof of 
Lemma 14.61 Using the Leibniz rule for difference quotients we obtain 

Recall that d s p ro — — d t p ro : so adding and subtracting gives 

If s = -S^AtlpWrdtPro + AtD^^rds Pr ■ 

Hence, 

\u A - u\ (D^ip + Lp s ) dX 

\ua - u\ S'i At '0w r (L»ip,. o - d t p ro )dX 

n| 

+ At \u A — it | D l _^\)ijj r d s p ra dX + \u A — u\ D_ipu) r p ro dX. 
J Ju 2 T J Jn% 

Finally, we use that 

D t _^ = 8+ t {t- V )-8l t {t-T). (6.6) 
Concerning the second term on the left-hand side, we apply p.8[) . The lemma now 
follows by sending e to zero, as in the proof of Lemma WM □ 



As seen by comparing Lemmas 16.21 and 15.31 there is one new term. To estimate 
this term we will use a result from [TJl p. 1853]. 



Lemma 6.3. Let u™ be the solution to (|6.1[) . Suppose the CFL condition (|6.2j) is 
satisfied. Then there exists a constant L such that 



AxJ2 \D+A(uf) - D+A(u%) \ < Ly/(m - n)At, for all m > 



n . 



j 



Estimate 6.1. Suppose (|6.2p is satisfied. Then 

At J J £)+sign E (A(u A ) - A(u))D-D+A(u A )tp dX 



Ax r V r o 
Proof. Integration by parts for difference quotients gives 

At JJ ^ Z^sign e (A(u A ) - A(u))D_D + A(u A )(pdX 

= -At j j sign E (A(u A ) - A(u))Sl At D-D + A(uA)D t _<pdX 



At / / sign e {A(u A ) - A{u))D'_D-D + A{u A )<pdX =: T x +T 2 , 



'li 2 

where we have used that 

Dl(D_D + i4(« A )v) = SlAtD-D+Aiu^D^tp + D'_D-D + A(u A )ip. 
Let us consider T± first. By the Leibniz rule for difference quotients, 

Ti = At JJ ^ sign £ (A(u A ) - A{u))S t _ At D_D + A{u A )S t _ M i>UrD t _Pr dX 
+ At JJ ^ sign e (A(w A ) - A(v))St_ M D-D + A(uA)D t _'il>WrPr dX 
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■ I . A I I . ,. 

Using equation (|6.6[) . 

|7i, 2 |<Af f |fil At D_D + A(u A )| (\8l t (t-v)\ + \SX t (t-T)\) dxdt<CAt, 

as \\D-D+A(uA(-,t))\\ L i^ is bounded independent of A and t ( [T3J Lemma 3.4]). 
Now, as in Lemma [4.111 



D-Pro < 



IP 



-l{|i-s|<r +At}(*) s ), 



and therefore 



\T t d < CAf ' + 9 A ^ / |5l A ,D_£) + yl(u A )| dxdt < C— (l + — 
r& Ju T r \ r Q 

Next, we consider T 2 . By Lemma 16.31 

|T 2 | < At [ f \D*lD-D + A(ua)\ dxdt 

Jv JR 

At r 

~ 2 Ai X ll^-^(«A(-,t))|| L1(R) ^ 
< 2TL 



Ax 



□ 



Proof of Theorem \6.1\ We start out from Lemma l6.2l with A(o~) = A(a)+rja, where 
A is the original degenerate diffusion function. By Estimate 16.11 

At JJ D\sign e (A(uA)-A(u))D-D + A(uA)<pdX 



Ax ro 



At 



CAt =: E 4 . 



Since all the estimates from Section 14. II apply, we obtain 

|ua - u| S At (t - r)uj r (x - y)pr {t - s) dX 



- JJ \ u &x~u\8 (t—v)w r {x — y)p ro {t-s)dX 

+ E 1 +E 2 + E 3 + E 4 , 

where Ei,E 2 , E 3 are defined respectively in (|4.9|l , (|4. 10[) , and (|5.5|l . Let us make 
the assumption that v = t m and r — t n for some m,n G N. Then the above 
inequality takes the form 



where 



K(t) :-- 



K(t n ) < re(t m ) + Ei+E 2 + E 3 + Ei, 



\v,A{x,t) - u(y,s)\w r (x - y)p ro (t - s)dydsdx. 



Applying Lemmas 14.121 and 15.11 and following the reasoning given in the semi- 
discrete case, we deduce 

\\u A (-,t n ) - u(-,t n )\\ L1(m) 

< \\u° A -u°\\ L i (m) + {L c + L d )t m 



40 



K. H. KARLSEN, N. H. RISEBRO, AND E. B. STORR0STEN 



2 (L c r + \u°\ Bvm r) + C(l + r + Axf (l + ^) 

Ax At / Ai\ „y/M 

C + C — 1 + r Q + — + C- 

ro r a V r J 



Ax\ 3 Ax 




< \\u° A -u° 

where Ld is the constant in Lemma 15.11 and L c is the constant from Lemma 14.121 
Let r = tq, Ax = 7- 3 and At = r 8 . It follows that 



\\u A (;t n ) -u(-,t n )\\ Ll{R} < \\u a A -u°\\ Ll . +CAx' 



~ 1/3 . 

Finally, we send r\ — > to conclude the proof of the theorem. □ 

7. Concluding remarks 

The added complexity of convection-diffusion equations versus conservation laws 
[2"2"] arises as a result of the need to work with an explicit form of the parabolic 
dissipation term. This is reflected in the fact that the rate of convergence is lowered 
to 1/3 (from 1/2 for conservation laws) due to Estimate 14.41 and Estimate 14.51 
The optimality of the ^ rate is an open problem. Concerning Section [6] (explicit 
schemes) , one may wonder if it is possible to remove the strengthened CFL condition 
At ~ Ax 8 ^ 3 (the usual one demands At ~ Ax 2 ). The difficulty is that the parabolic 
dissipation term is needed to balance the temporal error contribution as well as to 
carry out the doubling-of-the-variables argument, and this forces us to impose a 
stronger relation between At and Ax in order to appropriately control the temporal 
error contribution. We do not know if the condition At ~ Ax 8 / 3 is genuinely 
needed or is simply an artifact of our method of proof. Finally, we arc currently 
investigating the multidimensional case. For the semi-discrete scheme the main 
challenge seems to be the adaptation of Estimate 14.51 or more precisely to produce 
a multidimensional analogue of (|4.7|) . As an additional difficulty, Lemma RPl is not 
available in several space dimensions, see |12j . At the moment our multidimensional 
convergence rates are lower than in the one-dimensional case. 
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